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We examine the ground-state phase diagram of the t-J model in one dimension by means of 
the Density Matrix Renormalization Group. This model is characterized by a rich phase diagram 
as a function of the exchange interaction J and the density n, displaying Luttinger-liquid (LL) 
behavior both of repulsive and attractive (i.e. superconducting) nature, a spin-gap phase, and phase- 
separation. The phase boundaries separating the repulsive from the attractive LL phase as J is 
increased, and also the boundaries of the spin-gap region at low densities, and phase-separation at 
even larger J, are determined on the basis of correlation functions and energy-gaps. In particular, 
we shed light on a contradiction between variational and renormalization-group (RG) results about 
the extent of the spin-gap phase, that results larger than the variational but smaller than the RG 
one. Furthermore, we show that the spin gap can reach a sizable value (~ O.li) at low enough 
filling, such that preformed pairs should be observable at temperatures below these energy scales. 
No evidence for a phase with clustering of more than two particles is found on approaching phase 
separation. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm,74.20.Mn 



I. INTRODUCTION 

The t-J model constitutes together with the Hub- 
bard model a paradigm for the theoretical description 
of high temperature superconductors (HTS) since its 
derivation by Zhang and Rice 1 from a three-band Hub- 
bard (spin-fermion) model describing the copper-oxide 
planes present in HTS. It can also be derived in second 
order perturbation theory around U — oo, where U is the 
strength of the interaction in the Hubbard model 2 . Rs 
Hamiltonian is 

+J ■ - \ n i n i) > (!) 

where the operator fj (/j, CT ) creates (destroys) a fermion 
with spin a =1\ 4- on th° sit- They are not canonical 
fermionic operators since they act on a restricted Hilbcrt 
space without double occupancy. Si = fj a a a pf i ^ is the 

spin operator and i%i = f} a f^ a is the density operator. 
In all expressions a summation over repeated indices is 
understood. Furthermore, denotes nearest neighbor 
bonds. 

While the main interest on this model resides on its 
two-dimensional (2D) realization, it presents already in 
one dimension (ID) a number of very interesting features. 
In contrast to the ID repulsive Hubbard model, where 
only a Luttinger liquid (LL) phase for density n =/= 1 
and an insulating phase for n — 1 are present, the t- 
J model possesses a rich phase diagram, as shown first 
by M. Ogata et al. 3 . Interestingly, the phases display a 



correspondence to the ones present in HTS, like super- 
conductivity, spin gap and phase separation, albeit for 
values of J ft outside the range pertaining to HTS. How- 
ever, since unbiased results for the 2D model, that up to 
now could only be obtained by exact diagonalization 4 or 
by density matrix renormalization group (DMRG) , did 
not yet conclusively shed light on the different phases 
in the thermodynamic limit, the ID version presents a 
possibility of gaining insight into exotic phases like the 
spin-gap region. 

A further motivation for achieving an accurate de- 
termination of the phase diagram of this model is the 
possibility of realizing it with ultra-cold fermionic quan- 
tum gases 6 ' 7 . For parameters aiming at an emulation of 
HTS, still cooling techniques have to be implemented to 
reach the relevant energy scales (~ J/10, with J ~ 0.3i). 
On the other hand, as discussed in the following, the 
spin-gap phase appears here at much larger values of J 
(J ~ 2.5t), and hence this non-trivial phase is certainly 
much more accessible for experiments with degenerate 
quantum gases. 

The ID t-J model has been solved exactly only for 
J/t — > 0, where it is equivalent to the U — > oo Hubbard 
model 8-10 , and at the supersymmetric point J = 2t 11,12 . 
In both cases the model behaves as a LL 13 ~ 17 . For very 
large J/t the attractive interaction dominates against the 
kinetic energy and the system phase separates into hole 
rich regions and antiferromagnetic islands. Although 
the first phase diagram appeared almost twenty years 
ago, there are still issues to be clarified like the bound- 
aries of the spin-gap phase. In previous studies the 
existence of a spin gap has been deduced from exact 
diagonalization (ED) of small systems 3 , variational 18 , 
or projection 19 methods, transfer-matrix renormalization 
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group (TMRG), 20 and by a combination of renormaliza- 
tion group (RG) arguments with ED 21 . In particular the 
last study found a much larger region in density for the 
spin-gap phase. 

A main goal of the present work is to achieve a precise 
determination of the phase diagram, by performing finite- 
size extrapolations to the thermodynamic limit for the 
correlation functions and energy differences relevant for 
the different phases. The results are obtained using the 
DMRG method 22-24 , on lattices with up to L = 200 sites. 
We have to choose values of n and L such that the total 
particle number N = nL is an integer number. In most 
of the results we extrapolate to the thermodynamic limit 
using system sizes L = 40, 80, 120, 160, and 200. These 
values of L allow us a discretization in density of An = 
0.05 which is consistent with TV being an even integer, 
and hence, the ground-state corresponds to Sf ot = 0. 
All the results were obtained using at least 200 DMRG 
vectors, 4 sweeps, and a discarded weight of 10~ 8 . This 
translates into errors in energy of the order of 10 -8 and 
in correlation functions of the order of 10~ 5 at the largest 
distances. 
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FIG. 1. (color online). Phase diagram of the ID t-J model 
from DMRG for densities 0.1 < n < 0.9 and in the range 
< J < 4, where we set t = 1. n = N/L is the electronic 
density (TV is the total number of particles and L the number 
of lattice sites). Four phases are present: a metallic phase 
(M) or repulsive Luttinger liquid, a gapless superconducting 
(SC) phase, a singlet-superconducting phase with spin gap 
(SG + SS), and phase-separation (PS). The number given to 
each line stands for the value of the Luttinger parameter K p . 



II. PHASE DIAGRAM 

The phase diagram for the ID t-J model obtained by 
Ogata et al. 3 is based on ED on systems with up to 
16 lattice sites. They found three phases: repulsive LL 
phase (metal), attractive LL phase (superconductor) and 
phase-separation. At that time they suspected the exis- 
tence of a spin gap at low density but they could not 
prove its existence due to the limitation to small system 
sizes. Posterior works 18,19 found evidences of the spin 
gap using variational methods. However, Nakamura et 
al. 21 found based on an RG analysis that the spin-gap 
region is larger than expected. Here we present results 
from a direct measurement of the spin gap and an ex- 
trapolation to the thermodynamic limit. Details for it 
will be discussed in Sec. II B. 

Our results can be summarized in the phase diagram 
shown in Fig. 1. We obtained four phases: a metallic 
phase (M) or repulsive LL, a gapless superconducting 
(SC) region, a singlet-superconducting phase with spin 
gap (SG + SS), and phase-separation (PS), where the 
system separates into a hole-rich and an electron-rich 
part. The number given to each line stands for the value 
of the Luttinger parameter K p , with K p < 1.0 in M and 
K p > 1.0 in both superconducting phases. The determi- 
nation of K p will be discussed in the subsequent sections. 

In order to characterize the phases and to find the 
boundaries between them, we calculated directly the en- 
ergy gap to triplet excitations and measure the density- 
density correlation functions, 

Nij = (ni-rij) - (rii) (rij) , (2) 

the spin-spin correlation function, 

Sij — (SfSj), (3) 



the pairing correlation function 

Pii = (AjAj), (4) 

where 

for singlet pairing, and 

Al = 4 t 4i,t (6) 

for triplet pairing. Finally, we also considered the one- 
particle Green's function 

G1 3 = (/i/j). (7) 

The corresponding structure factors are obtained by 
Fourier transformation, 

1 L 

X ^ = ^E ^~ X]) Xr 3 , (8) 

Although the systems considered lack translational in- 
variance due to open boundary conditions, for the large 
system sizes considered here, we could not observe any 
artifact introduced by this procedure. 

A. Metallic phase 

In order to characterize this phase, we compute the 
Luttinger parameter K p , with K p < 1 (K p > 1) for a 
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repulsive (attractive) interaction, and K p = 1 for the free 
case. In oder to obtain K p , we consider the limit k — > 0, 
where the structure factor for the density correlations 
displays a linear behavior with a slope proportional to 

TS 13,25,26 

1v p > 



N(k) -> K p \k\a/n for k -> 0, 



(9) 



that results from Fourier transforming the first term in 
Eq. (12) below. Here a is the lattice constant (we set 
a = 1). Figure 2 shows N(k) for n — 0.5, J = 2.0 and 
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K p as function of J for different densities n is 
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We repeated this procedure for different values of n 
and J 

plotted in Fig. 4. Note that K p 
for all densities, which is in agreement with the results 
obtained for the U/t — > oo Hubbard model 8 . It can also 
be observed that, for K p > 1, K p increases quite fast 
with the interaction constant J and it actually should 
diverge in the phase-separated region. 




FIG. 4. (color online). K p as function of J for different den- 
sities n. 



FIG. 2. (color online). Structure factor N(k) of the density- 
density correlation function for n = 0.5, J = 2.0 and L = 40, 
80, 120, 160, and 200. 

L = 40, 80, 120, 160, and 200. We observe a clear linear 
behavior for small k, with N(k = 0) = due to the 
conservation of the total particle number. 

Although N(k) appears to be almost independent of 
the lattice size, a more precise value of the slope is ob- 
tained by extrapolating the value of N(k) at the point 
k = 27r/40, that is the smallest wavevector in our small- 
est system, to L — » oo. Using this last value and 
N(k = 0) = we obtain the slope and then we can ex- 
tract K p in the thermodynamic limit using Eq. (9). This 
extrapolation is shown in Fig. 3. 
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FIG. 3. (color online). Extrapolation to the thermodynamic 
limit of N(k = 2tt/40) for n = 0.5 and J = 2.0. 



The critical exponents K p at the supersymmetric point 
,1 = 1 and for all densities were exactly obtained by 
means of the Bethe ansatz 27 . In Fig. 5 we compare our 
DMRG results with this exact solution, and observe a 
very good agreement between both. The deviations at 




FIG. 5. K p as function of the density n for J = 2 (supersym- 
metric point). 

very low (very high) densities point to the necessity of 
having larger systems for such very dilute cases, with a 
rather small number of particles (holes) to be able to 
properly describe a phase. They lead however, to barely 
noticeable shifts in the phase diagram. 

From the data set presented in Fig. 4 we can extract 
all the points which fulfill K p (n,J) = const. These are 
curves which separate regions with different Luttinger 
parameters. The different regions and curves are plotted 
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in Fig. 6. The red (dashed) line (K p (n,J) = 1) denotes 
the boundary between the metallic phase and the super- 
conducting phases. This line and few others were also 
plotted in the phase diagram (Fig. 1). Note that the 
density of lines increases with J showing the fast growth 
of K p . 
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FIG. 7. (color online). Energy comparison between systems 
containing four, six and eight particles and a gas of singlet 
bound pairs. We observe a region 2.0 < J < 3.0 where all 
energies are the same within an error of 10~ 5 . This opens 
the possibility of the formation of a gas of singlet bound pairs 
and consequently of the existence of a spin gap at very low 
densities. 



FIG. 6. (color online). Regions with different Luttinger pa- 
rameters K p (n,J). Each curve corresponds to a constant 
value of K p . The red (dashed) line (K p (n,J) — 1) denotes 
the boundary between the metallic phase and the supercon- 
ducting phases. Note that the density of lines increases with 
J showing the fast growth of K p . 



B. Singlet-superconductivity and spin-gap phase 

The possibility of a region with a spin gap was first an- 
alyzed by Ogata et al. 3 for the low density limit, where 
a gas of singlet bound pairs may form. They compared 
the ground state energy of a system containing four par- 
ticles to twice the energy of a system with only two par- 
ticles. The energy for the last situation is 2(— J — 4/J), 
where the expression — J — 4/J is obtained by solving 
exactly the Hamiltonian (1) for two particles. We com- 
pare here the ground-state energy per particle for 4, 6, 
and 8 particles obtained numerically with DMRG to that 
of a gas of bound pairs (Fig. 7), using 1000 lattice sites. 
We observe a region, 2.0 < J < 3.0, where energies are 
the same within an error of 10~ 5 . At least from energy 
considerations, this is an indication of the possibility of 
the formation of a gas of singlet bound pairs and conse- 
quently of the existence of a spin gap at very low den- 
sities. Moreover, no evidence for the formation of more 
complex entities can be seen. 

A more precise estimate can be obtained by measuring 
directly the spin gap AE$- The spin excitation energy 
from a singlet to a triplet state is given by the energy 
difference 

AE S = E (N, S? ot - 1) - E (N, S z tot = 0), (10) 

where the subindex means that we take the lowest en- 
ergy level with given quantum numbers ./V and S% ot . In 



order to go to the thermodynamic limit we plot AE$ as 
a function of 1/L and we extrapolate to 1/L using 
L = 40,80,120,160, and 200. Figure 8 shows AE S vs 
1/L at n = 0.2 and for J = 2.0, 2.7 and 2.9. The extrap- 
olations to the thermodynamic limit are performed with 
third-order polynomials. While for J = 2, where the sys- 
tem is still in the metallic regime the gap extrapolates to 
zero, for J = 2.7 and 2.9 a finite gap can be clearly seen. 
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FIG. 8. (color online). AEs vs 1/L for n = 0.2 and various 
values of J. For J = 2, where the system is still in the metallic 
phase, the spin gap extrapolates to zero in the thermodynamic 
limit. 

Proceeding in the same manner for different values of n 
and </, we can obtain the spin gap in the thermodynamic 
limit (Fig. 9). For J < 2 (metallic phase) we observe 
that the spin gap is zero for all densities. For J > 2 
a finite spin gap emerges that increases as the density 
diminishes. For definiteness, J c is defined in our case 
as the value of J for which AEg > 10~ 4 t, this value 
being the the range on which AE$ fluctuates around zero 
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before it definitively increases as a function of J for a 
given density. In this manner we have obtained the lowest 
boundary of the spin-gap phase in the phase diagram 
(Fig. 1). 
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FIG. 9. (color online). Spin gap AE$ in the thermodynamic 
limit as a function of J and for different densities n. For 
J < 2 (metallic phase) the spin gap is zero for all densities. 
For J > 2 a finite spin gap emerges with a value that increases 
for diminishing densities. We can observe a small but finite 
spin gap up to n = 0.55. 



We present also in Fig. 10 the spin gap as a function 
of n for J = 2.1 — 2.8. Note that AE$ smoothly closes 
to zero when the density is increased. AE$ attains its 
largest values as J increases for the limit of vanishing 
densities, reaching AE$ ~ J/20. Since in one dimension 
superconductivity can only set in at temperature T = 0, 
such a finite value of the gap gives the energy scale at 
which pairs form, signaling the existence of preformed 
pairs in this regime. A further increase of J leads to 
phase separation, that we discuss next. A detailed anal- 
ysis of the transition to phase separation, centering on 
the possibility of clusters with more than one pair is pre- 
sented in Sec. HID. 
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FIG. 10. (color online). Spin gap AE$ in the thermodynamic 
limit as a function of n for J = 2.1 — 2.8. 



C. Phase separation 



In this phase the attraction among the particles is 
so strong that they start to form antifcrromagnetic do- 
mains, such that the system separates into particle- and 
hole rich regions. In the limit J — > oo all the particles 
join in a single island, which can be described by the 
Heisenberg model forming an electron solid phase, a de- 
nomination proposed by Chen and Moukouri 28 , where 
the kinetic fluctuations are strongly quenched and only 
spin fluctuations remain. We first consider the inverse of 
the compressibility that vanishes at the onset of phase- 
separation. 
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FIG. 11. (color online). Inverse of the compressibility k 1 as 
a function of J for n = 0.1 — 0.9 with An = 0.05. 



At zero temperature the expression for the inverse com- 
pressibility is given by 



-1/ \ 2 <9 2 eo(n) 
k [n) = n 



dn 2 

>(n + An) + e(n - An) - 2e(n)] 
An 2 



, (11) 



where eo(n) = Eq/L is the energy density per site, and 
the second line gives the approximation for finite (An = 

0. 05) changes in the density. 

For the extrapolation of e(n) we use L = 
40, 80, 120, 160 and a third-order polynomial fitting. Fig- 
ure 11 shows k _1 vs J for different densities. At low 
densities k _1 is vanishingly small, making the extraction 
of J c very difficult in that region. In order to see more 
clearly the critical value J c {n) where k _1 vanishes, we 
display a zoom of Fig. 11 in Fig. 12. In this manner we 
found the boundary of the phase-separated phase in Fig. 

1. Note that, in comparison to other studies 3 ' 18,19 , the 
phase separation boundary is shifted to higher values of 
J. 
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FIG. 12. (color online). Zoom of Fig. 11. The points (n, J c ) 
where k" 1 — define the boundary of the phase-separated 
phase (infinite compressibility). 



III. CORRELATION FUNCTIONS AND 
DENSITY IN REAL SPACE 

A. Structure factors 

In order to provide a more detailed characterization 
of the different phases in the phase diagram determined 
in the previous section, we consider here the structure 
factors Eqs. (2) - (8). They are shown in Figs. 13-16 for 
L = 200 and for different values of n and J. 

Figure 13 displays the structure factor N(k) for the 
density-density correlation function, that for a LL is as 
follows 8 : 



(n(r)n(O)) = 



cos{2k F r) _ 3/2 

A l ,!+*„ 111 ( r ) 



+A2 cos(4fc^r)r 



(12) 



where A± and Ai depend in general on the parameters 
of the model. For each density we have chosen values of 
J within one of the phases. For J = 1, a LL phase is 
realized for all densities, such that a 4kp (kp = mr/2) 
anomaly typical for a repulsive LL can be observed. As 
shown by Eq. (12), it is strongest for the smallest val- 
ues of K p , that are essentially achieved for all densities 
for J < 1, as shown by Fig. 4. As J increases, the 4kp 
anomaly is suppressed, and a 2kp cusp is formed, signal- 
ing a 2kp charge density wave due to the tendency to- 
wards pairing that is enhanced by J 13 ' 29 . However, going 
closer to the boundary to phase separation (see the curves 
for J = 3 in Fig. 13 (b), and for J = 3.15 in Fig. 13 (c)), 
the singularity is rounded, and increasing J even more, 
such that phase-separation is reached, leads to the devel- 
opment of a singularity in N(k) around k — 0. That is, 
the system starts to develop an instability towards long- 
wavelength charge fluctuations signaling the appearance 
of phase-separation. However, since in the finite system 
simulations, N(k = 0) = for all cases due to the con- 
servation of total particle number, such a singularity can 
only be followed up to the smallest non- vanishing value of 
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FIG. 13. (color online). Structure factor N(k) for the density- 
density correlation function for L = 200 and for different val- 
ues of n and J. The location of 2/cf and 4fci? correspond to 
the ones resulting from folding them back to the sector k > 
in the first Brillouin zone. 



momentum for a given system size. Except for the cases, 
where the system enters the phase separated phase, N(k) 
goes continuously to zero, as k — ¥ 0, such that K p can be 
extracted, as discussed in Sec. II. In both superconduct- 
ing phases, K p > 1, as expected, and increases as one 
goes deeper into the spin-gap phase. 

In Fig. 14 the structure factor S(k) for the spin-spin 
correlation function, is shown. The spin-spin correlation 
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FIG. 14. (color online). Structure factor S(k) for the spin- 
spin correlation function for L — 200 and for different values 
of n and J. 



function is given by 8 



(S z (r)S z (0)) 



(7rr) 2 



D 



1 CC ^ln^ (r ).(13) 



The tendency to antiferromagnetism in the t-J model can 
be observed at J = 1 for all densities, as revealed by a 
sharp cusp at 2kp in S(k), corresponding to quasi-long 
range order in the magnetic channel. However, on in- 
creasing J, such that the system enters the supercon- 
ducting phase, the sharp peak is suppressed. For den- 
sities below n = 0.6, a spin gap develops, such that on 
entering this Luther-Emery (LE) phase, due to the ex- 
ponential decay of the spin-spin correlation function, the 



singularity at 2kp is completely suppressed 13 . In this 
case, the corresponding structure factor has a quadratic 
behavior at small fc's 18 . This fact can be clearly observed 
only at J = 2.5 and n — 0.1, where the spin gap is well 
developed. As shown in Fig. 9, the spin gap at J = 3 for 
a density n — 0.5 is extremely small (AEs ~ 10~ 3 ), such 
that extremely large systems would be required to show 
such a behavior. On increasing J a maximum at k — 7r 
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FIG. 15. (color online). Structure factor Ps(k) for the singlet 
pair-pair correlation function for L = 200 and for different 
values of n and J. 

develops due to the tendency to form antiferromagnetic 
islands in the phase-separated region. However, while 
the charge structure factor already shows clear signals of 
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phase-separation (Fig. 13), S(k) displays only a broad 
maximum around k = ir not yet indicative of antiferro- 
magnetic quasi-long range order. A larger value of J is 
necessary in order to achieve such a state, as suggested 
by the curves corresponding to the largest value of J in 
Fig. 14. Hence, the onset of phase-separation and the 
formation of an antiferromagnetic island do not occur si- 
multaneously, as already observed in an earlier quantum 
Monte Carlo study 29 . We will discuss the formation of 
antiferromagnetic islands in more detail in Sec. Ill D. 

In Fig. 15 the structure factor Ps(k) for the singlet 
pair-pair correlation function is shown. The correspond- 
ing correlation function in the LL sector is given by 13 

(A t s (r)A s (0)) = C r-( 1+1 /^>) 

+Ci cos(2k F r)r- {K » +1/K »\ 
and in the LE phase by 

(At (r)A s (0)) = C' r^' K ? 

+C[ cos(2k F r)r~ 



-{K p + l/K„) 



(14) 



(15) 



where we have ignored logarithmic corrections. Noticing 
the different powers appearing in Eqs. (14) and (15), and 
taking into account that K p > 1 in the superconduct- 
ing phases, it can be seen that while in the LL sector 
Ps(k) does not present a divergence, it will have a di- 
verging contribution at k — in the LE case. Such a 
behavior can be seen in Fig. 15. Within the spin-gap 
phase (Fig. 15 (a), J = 2.5 and Fig. 15 (b), J = 3), 
a strong increase is observed for k — > 0, indicating the 
onset of quasi-long range order in this channel, (singlet- 
superconductivity) . While increasing J and n also an en- 
hancement of Ps (k = 7r) is produced in the LL phase, the 
curves display only a rounded maximum in this case. We 
observe also &2kp anomaly in Pg(fc) that correlates with 
the 2kp singularities in N(k), displayed in Figs. 13. We 
will discuss these correlation functions further in Sees. 
Ill B and III C, where they will be confronted with other 
possible orderings like triplet pairing and charge-density 
wave (CDW) formation. 

In Fig. 1G the momentum distribution function n(k) 
is shown. In the metallic phase an edge develops in 
n(k) that evolves into a singularity in the thermodynamic 
limit, and hence, defines a Fermi surface, in the way ex- 
pected for a LL. A weak anomaly at 3k F can be observed, 
as revealed by the derivative of n(k) shown in the insets in 
Fig. 16. As was analytically shown by Ogata and Shiba 30 
for the Hubbard model, this singularity in n(k) is related 
to the fact that one electron close to kp can be excited to 
a state close to 3kp together with an electron-hole pair 
excitation having its momentum near ~2kp. In Fig. 16 
the 3kp anomaly disappears as J increases, in agreement 
with previous DMRG studies 28 . 

For n = 0.1 and J = 2.5, which corresponds to a point 
in the spin-gap phase, we observe that the Fermi sur- 
face is destroyed. For the higher densities in Figs. 16 
(b) and (c), such a flattening of n(k) around kp can- 
not be observed because the spin gap for such parame- 
ter values is very small. An interesting feature seen for 
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FIG. 16. (color online). Momentum distribution function 
n(k) for L = 200 and for different values of n and J. 



such densities is the increase of n(k) for k > kp, con- 
trary to what is expected in a weakly correlated metal. 
Such a behavior is indicative of the presence of spectral 
weight below the Fermi energy for wavevectors k > kp, as 
was already observed in previous quantum Monte Carlo 
simulations 31 . On the contrary, for n = 0.1 and J = 2.0 
(where K p « 1.0) n(k) has a shape closer to the momen- 
tum distribution in the free case. 
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B. Gapless superconducting phase 

In this section we examine closer the pairing correla- 
tion functions both for singlet as well as for triplet pair- 
ing. The long-distance behavior of the correlation func- 
tion for singlets was already given in Eq. (14) and for 
triplet pairing is as follows 13 ' 32 : 

(A T (r)A T (0)) ~ A)^ (1+1/A ' p) 

+£>! cos(2k F r)r- iK » +1 / K » +2 \ (16) 

where the operator At corresponds to triplet supercon- 
ducting pairing (see Eq. (6)). Here we have ignored log- 
arithmic corrections. 

In the gapless superconducting phase, since there is no 
gap to triplet states, and the leading power-law is the 
same for both singlet and triplet channels, the question 
may arise about the relative strength of singlet- (SS) and 
triplet-superconductivity (TS), as already pointed out by 
Pruschke and Shiba 32 . 




X 

FIG. 17. (color online). Singlet Ps(x) and triplet Pt{x) cor- 
relation functions for n — 0.5, J = 2.6 and L — 200. For the 
singlet case we plot |Ps(x)|. The straight lines correspond to 
power-laws determined by K p . Both the singlet and triplet 
channel have the same exponent. 

Figure 17 shows the behavior of both correlation func- 
tions in real space as compared to the asymptotic forms 
given in Eqs. (14) and (16). In Fig. 17 (a) a compari- 
son is given with Eq. (14), where the constants Co and 



C\ were adjusted through a least-square fit, while the 
value for K p was taken from the determination detailed 
in Sec. II A (K p = 1.12 for the parameter values in Fig. 
17). In order to be able to display the power-law be- 
havior of P$(x), we use a doubly logarithmic scale, and 
actually plot the modulus of Ps(x), since the 2Uf oscil- 
lations lead to sign changes of the correlation function. 
The dashed line through the maxima makes furthermore 
evident, that the power-law decay can be well described 
by the power of the first term in Eq. (14). Deviations 
are observed for x > 60, possibly due to boundary ef- 
fects. On the other hand, Fig. 17 shows that the decay 
of Pt can be well described by the same power-law as 
for Ps(x), with deviations that start at x > 30. While 
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X 



FIG. 18. (color online). Singlet Ps(x) and triplet Pt(x) 
correlation functions for n = 0.5, J = 2.6, L = 200 and 
m = 150, 200, 250. We observe at long distances that on a 
double-logarithmic scale, Pt(x) becomes more linear upon 
increasing m, i.e., the spurious exponential decay introduced 
by the DMRG cutoff decreases. Note, however, that Ps{x) 
becomes less linear at long distances when increasing m. We 
associate this behavior to the open boundary conditions used. 



there are certainly boundary effects, Fig. 18 shows that 
their incidence on the correlation functions changes by 
increasing the precision of the DMRG runs. There we 
display the results for both correlation functions when 
the number m of states kept in the reduced density ma- 
trix is increased. While Pt(x) approaches the power- 
law at larger distances, Ps(x) departs from it. At the 
highest accuracy used, both depart from the predicted 
power-laws for x > 50. 

In Fig. 19 we also show both singlet and triplet pair 
correlations in momentum space. The amplitude of 
Ps(k — 0) for the singlet case is one order of magni- 
tude larger than Pr{k = 0) for the triplet case. Since 
Ps(k = 0) gives the number of pairs with momentum 
zero, it is clear that singlet superconductivity dominates. 
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FIG. 19. (color online) . Triplet and singlet pair-pair structure 
factors for n — 0.5, J = 2.6 and L — 200. 



C. Competing orders in the spin-gap phase. 
Charge density wave vs singlet superconductivity 

The peaks of Fig. 13 (a) at 2k F and Fig. 15 (a) at k = 
for J = 2.5 show that CDW and SS are competing orders 
in the spin-gap phase. On entering this phase, electrons 
pair into singlets, such that the spin-spin and triplet pair- 
pair correlation functions are exponentially suppressed. 
While the singlet pair-pair correlation function has the 
long-distance behavior given by Eq. (15), the density- 
density correlation function in the LE phase behaves as 13 



(n(r)n(O)) = A r 2 + A x cos(2k F r)r 



-K 



(17) 



Since in the spin-gap region K p > 1, SS is the dominant 
order. The difference between both correlation functions 



i(r 
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FIG. 20. (color online). Maxima of Ps(x) and —N(x) in real 
space for L = 200, n = 0.2 and J = 2.5 on a log-log scale, 
the power-laws are determined by K p through Eqs. (15) and 
(17). 

is displayed in Fig. 20, where we consider them at n = 0.2 
and J = 2.5, i.e. deep in the spin-gap region. Due to the 
2kp oscillations, only the maxima of the functions are 
plotted. The different lines (dashed one for Ps(x) and 
full and dashed-dot lines for N(x)) correspond to the 
powers appearing in the first term of Eq. (15), and the 



powers appearing in Eq. (17) for N(x). The latter shows 
a crossover from a behavior at short distances dominated 
by the first term in Eq. (17) to the long-distance behav- 
ior determined by the power of the second term. For the 
parameters considered here we obtained from N(k — > 0), 
as given by Eq. (9), K p = 1.48. In this case Ps(k = 0) 
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FIG. 21. (color online). Scaling of Ps(k = 0) in the thermo- 
dynamic limit for n = 0.2 and J = 2.5 (spin-gap phase). 

diverges in the thermodynamic limit, as shown by Fig. 21 
with an exponent determined by K p , however, the sys- 
tem being one-dimensional, no true long-range order is 
present, since Ps(k = 0)/L vanishes in the thermody- 
namic limit. Nevertheless, as shown in Fig. 10, the spin 
gap increases as the density decreases, reaching a sizable 
value on approaching phase separation (AE$ ~ 0.15t). 
Hence, singlets bind at a finite temperature scale below 
the spin gap. 

Another interesting aspect of this phase is that due 
to open boundary conditions, a change in the periodic- 
ity of the Friedel oscillations in the density can be di- 
rectly observed, as shown in Fig. 22. This can be seen 
as a 2kp CDW coexisting with a superconducting state 
in the spin-gap phase. However, following the arguments 
related to Fig. 7, and the fact that the leading singu- 
larity is related to SS, these density oscillations can be 
understood as due to bound pairs, which because of the 
constraint of the t-J model behave as hard-core bosons. 
Hence, the density of the hard-core bosons in one dimen- 
sion will show the same oscillations as the density of the 
equivalent Jordan- Wigner fermions, corresponding to the 
number of pairs in the system. As an example we show 
in Fig. 22 the density profiles for a density n — 0.1 and 
values of J corresponding to the repulsive LL (J = 1) 
and the spin-gap (J = 2.5) phases, and two values of J 
on entering the phase separation, for two different sys- 
tem sizes (Fig. 22 (a): L = 80 and Fig. 22 (b): L = 160). 
In both cases, the number of oscillations is halved, cor- 
responding to Friedel oscillations for half the number of 
particles. 

Since the coupling driving the pairing mechanism is the 
nearest neighbor exchange J, it could be argued 3 that be- 
fore entering phase separation, clustering in groups with 



11 



<W=1.0 

J=2.9 
J=3.0 




FIG. 22. (color online). Density in real space p(x) for n — 0.1 
and (a): L = 80, (b): L — 160. In the metallic regime (J = 
1.0) we have one peak per particle due to Friedel oscillations. 
When J is increased the particles start to form pairs. This 
happens in the spin-gap region (J = 2.5). Increasing J even 
more the particles tend to be confined in a region smaller than 
the system size (Phase separation). 



more than two fermions may exist. In order to answer 
this question we observe how the modulations of p(x), 
displayed in Fig. 23, change when we increase J by very 
small amounts before phase-separation is reached. There, 
we clearly see how a loosely bound cloud forms at the bor- 
der of the particle rich region, but pairs remain as such in 
the middle of that region. Starting with 16 electrons on 
160 sites, we can see the formation of pairs on entering 
the spin-gap region, such that 8 pairs are clearly visible 
(Fig. 23 (a)). Going up to the point where density oscil- 
lations corresponding to two pairs are still visible (Fig. 
23 (d)), we see that the rest of the pairs merged on the 
sides of the particle rich region, until phase-separation is 
reached. Note that the density oscillations in Figs. 23 are 
on a very small scale and are not noticeable at the scale 
used in Figs. 22. Once in the phase-separated region, a 
density distribution results, where a cloud at the center 
of the system appears, leaving an appreciable number of 
empty sites, as shown in Figs. 22. Hence, clusters with 
more than a pair do not form a uniform phase, and in 
particular, a state with four bound electrons does not 
form. Still inside the phase-separated region, apprecia- 
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FIG. 23. (color online). Density in real space p(x) for n = 0.1, 
L — 160 and different values of J. In Fig. 23(a) we observe 
the pairing of particles. In the other figures only particles at 
the boundaries of the particle-rich region merge into one loose 
cloud but the particles in the middle do not cluster further. 



ble changes of the density as a function of J take place. 
This will be the subject of the next section. 
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D. Phase separation and electron solid phase 
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FIG. 24. (color online) Magnetic structure factor S(k) in the 
phase-separated region for values of J close to the boundary 
and at the value, where a peak at k — ir emerges. Insets: the 
corresponding density profiles 

As shown in Fig. 22 at values of J just after the phase 
boundary (J = 3), particles merge into a single island, 
occupying only part of the available space. However, as 
shown by the magnetic structure factor in Fig. 14, for val- 
ues of J inside phase-separation, and close to the bound- 
ary, a broad maximum is seen at k — n, but no sharp 
peak indicative of the formation of a spin chain. There- 
fore, the naive picture of a compact region does not apply 



yet. Figure 24 shows S(k) at three different densities for 
the values of J where the system enters phase-separation 
and the one at which a sharp maximum can be observed 
at k = tt. The insets display also the corresponding den- 
sity profiles, making evident that when S(k) has a sharp 
maximum at the antiferromagnetic wavevector, an island 
with density n = 1 is formed. Unfortunately, since the 
formation of such islands implies that there are many 
almost degenerate states very close to the ground-state, 
namely those connected by translation, is is not possible 
to perform a careful finite-size analysis to determine the 
boundary to such a phase in the thermodynamic limit. 
In fact, in Figs. 24 it can be seen that the density profiles 
reaching n = 1 break spontaneously reflexion symmetry 
about the central bond, an artifact due to the many de- 
generate states mentioned above. 
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FIG. 25. (color online) Phase diagram of Fig. 1 including the 
electron solid phase (ES). The circles denote the position at 
which a peak in S(k) appears at k = n, and at the same time, 
the phase separated island reaches a density n = 1. Triangles 
reproduce the results of Ref. 21 and the dashed-dot line those 
of Ref. 19. 

The difference in density profiles shown in Figs. 24 
makes also evident that by having an expanded cloud, 
most probably with pairs inside, as discussed in the pre- 
vious section, still kinetic energy is present in the cloud 
that is quenched only when J reaches a large enough 
value. Since we are limited in the sizes of the systems 
we can simulate, we cannot determine with confidence a 
critical value of J, where the system forms an electron 
solid 28 . It was argued 28 that in this case J c should be 
independent of the band filling, n, since the only require- 
ment is to form an island with n = 1. We observe in 
our simulations a weak dependence on the density, but 
our results are limited to L — 80, the largest size where 
meaningful results can be obtained in this region. Plac- 
ing the location of the value of J, where the island with 
n = 1 forms for L = 80 at different densities, leads to an 
almost density independent value, as shown in Fig. 25. 
The largest deviations are observed at the lowest den- 
sities. Although the number of particles is possibly too 
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low to make a definite statement, it could be argued that 
at low densities the formation of the electron solid leads 
to an appreciable loss of kinetic energy, that has to be 
compensated by a larger value of J c . The same can be 
argued in the high density region, such that being limited 
to small systems may lead to lower values of J c , since, 
as displayed in Fig. 24 (c), for very few holes the loss 
in kinetic energy with open boundary conditions can be 
negligible. 

Figure 25 also displays the differences in the results 
for the onset of the spin-gap region that were obtained 
on the basis of a variational calculation combined with 
the power method 19 on the one hand, and a RG analysis 
complemented by exact diagonalizations of systems with 
up to 30 sites 21 , on the other hand. While the agreement 
with the latter results is in general quite good, the largest 
deviations appear as the onset line approaches the phase- 
separation line. This can be possibly due to the fact that 
the RG analysis did not take explicitly into account the 
appearance of phase-separation, that in general leads to 
a strong finite size dependence. On the other hand, a 
precise determination with DMRG becomes extremely 
demanding at low densities (n < 0.1), due to the need of 
very large systems in order to have a large enough number 
of particles to faithfully represent a phase. The large 
disagreement with the results of Ref. 19 could point to 
the difficulties of the power method to deal with fermionic 
systems at high density, due to the strong increase of the 
dimension of the required Hilbert space. 



IV. SUMMARY 

We have revisited the phase-diagram of the one- 
dimensional t-J model and determined, on the basis of 
finite-size extrapolations of results obtained with DMRG, 
the boundaries between the known four phases: metal 
(M), singlet-superconductivity with spin gap (SG+SS), 
gapless superconductivity (SC), and phase separation 
(PS) (see Fig. 1) . The most controversial issue was re- 
lated to the boundary between SC and SG+SS, where 
appreciable differences were present between the re- 
sults from variational methods 18 ' 19 , and results based 
on renormalization group 21 . The highest densities at 
which the spin-gap phase was predicted was n ~ 0.4 
and n ~ 0.8, respectively. In our case it corresponds 
to n ~ 0.6. The boundary between M and SC was de- 
termined by extracting the Luttinger liquid anomalous 
dimension K p from the slope of the structure factor for 
density correlations in the limit k — > 0, extrapolated to 
the thermodynamic limit. The extrapolations were per- 
formed using system sizes L = 40, 80, 120, 160, and 200. 
The opening of the spin gap was directly determined by 
examining the gap to the lowest triplet state, and again 
extrapolating to the thermodynamic limit. Finally, the 
boundary to PS was determined by extrapolating the in- 
verse compressibility. 

We further characterized the different phases through 



correlation functions for density, spin, singlet and triplet 
superconductivity, and their corresponding structure fac- 
tors. The correlation functions could be consistently de- 
scribed by the determined values of K p . We also consid- 
ered the momentum distribution function, that at high 
densities shows a peculiarity not possible in a single band 
conventional metal. Apart from detecting the 3kp singu- 
larity, at high enough densities, n(k) increases for k > kp. 
This corresponds to spectral weight for energies below 
the Fermi energy but k > kp, as observed in previous 
quantum Monte Carlo simulations 31 . 

In the gapless superconducting phase the Luttinger 
parameter K p is larger than one and both SS and TS 
correlation functions decay with the same critical expo- 
nent. The amplitude of the structure factor for singlet 
pairing clearly dominates over that for triplet pairing, as 
expected (Fig. 19). A very good description of the cor- 
relations in real space is obtained by the forms given by 
Luttinger liquid theory (Fig. 17), and with power-laws 
consistent with K p determined from the structure factor 
for the density-density correlation function. It describes 
both the decay of the correlation function for the singlet 
as well as the triplet channel. 

In the spin-gap phase with singlet-superconductivity 
we observed that the ground state energy compares very 
well to the ground state energy of a gas of singlet bound 
pairs in a region 2.0 < J < 3.0 at very low densities (Fig. 
7). Correspondingly, a,2kp singularity due to pairing can 
be seen in the structure factor N(k) (Fig. 13 (a) and (b)). 
In fact, the density profile of systems with open boundary 
conditions shows in this phase modulations that can be 
viewed as pairing in real space (Fig. 22). We would like 
to also remark that the energy scale of the spin gap can 
reach values ~ O.lt, and hence their formation should be 
experimentally accessible at finite temperatures. 

The expectation of clustering of electrons beyond pairs 
close to the boundary to phase separation as a possible 
phase is not supported by our calculations. While in- 
creasing J at low density n, we monitored the particle 
density modulations in real space (Fig. 23). We see that 
between the LL and the spin-gap phase a pairing of par- 
ticles in real space in fact occurs. However, on increasing 
J the particles build a loose cloud only at the boundaries 
between the particle-rich region and the hole region. In 
the center of the particle- rich region still modulation cor- 
responding to pairs can be seen. 

On entering phase-separation, the following features 
are present: infinite compressibility (Fig 12), Luttinger 
parameter K p — > oo (Fig. 4), and, confinement of parti- 
cles in real space to a region smaller than the available 
space (Fig. 22). In comparison to other studies 3,18,19 , 
the phase-separation boundary is pushed up to higher 
values of J. As observed in early quantum Monte Carlo 
simulations 29 we see that the onset of phase separation 
and the formation of a single antiferromagnetic Heisen- 
berg island do not occur simultaneously. Although a 
detailed finite-size scaling was precluded by metastabil- 
ity problems in the numerical implementation, we deter- 
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mined an approximate boundary to the appearance of 
such islands, termed previously electron solids 28 . Figure 
25, containing these data displays the region, where such 
islands occur. 
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